AJ in press 



High Redshift X-Ray Selected Quasars: 
CXOCY J125304.0-090737 joins the club 1 

Francisco J Castander 2 ' 3,4 ' 5 , Ezequiel Treister 2,3 , Thomas J Maccarone 6 , Paolo S Coppi 2 , 
Jose Maza 3 , Stephen E Zepf 7 and Rafael Guzman 8 

ABSTRACT 

We present a new X-ray selected high redshift quasar CXOCY J125304.0- 
090737 at z — 4.179, discovered by the Calan-Yale Deep Extragalactic Research 
(CYDER) Survey. This quasar is the fifth X-ray selected high redshift radio quiet 
quasar {z > 4) found so far. Here, we present its observed properties which are 
characterized by its relative optical and X-ray faintness, its X-ray hardness and 
its X-ray strength compared to optically selected quasars at high redshift. 

We also compare the X-ray selected high redshift radio quiet quasars to their 
optically selected counterparts. We find that the optical to X-ray spectral slope, 
a ox , is statistically harder (more X-ray luminous) for the X-ray selected radio 
quiet quasars than for the optically selected ones. This result, given the differ- 
ent range of rest frame ultraviolet luminosities studied and the selection of the 
samples, is consistent with the previously found correlation between X-ray and 
rest frame ultraviolet luminosities and would extend that result to a much wider 
luminosity range at high redshift. Finally, we discuss the prospects of unveiling 
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the quasar luminosity function at high redshifts using X-ray surveys. The dis- 
covery of a high redshift object in the first field of our survey program provides 
suggestive evidence that X-ray selected surveys may identify more such objects 
than would be expected from an extrapolation of the optical luminosity function. 

Subject headings: galaxies: active, quasars: individual, X-rays: galaxies 

1. Introduction 

The formation of the first structures in the Universe is one of the most interesting 
subjects in today's astronomy. Because quasars are intrinsically very luminous, they can 
be observed up to very high redshift and can therefore probe early epochs in the universe. 
Recent observations have found considerable numbers of quasars at high redshift (Anderson 
et al. 2001; Schneider et al. 2002), with some quasars being at redshifts ~ 6 (Fan et al. 
2001b), when the Universe was less than one Gyr old. 

The earliest quasar surveys were performed in the radio (e.g., Edge et al. 1959; Bennett 
1962; Schmidt 1963) and the optical (e.g., Sandage 1965). While radio surveys have continued 
(e.g., White et al. 2000), in recent years, greater attention has been placed on optical and X- 
ray surveys. Optical surveys are efficient finding quasars as they can search large areas of the 
sky in reasonable amounts of time. At high redshift (in this paper meaning redshifts z > 4), 
the Sloan Digital Sky Survey (SDSS; York et al. 2000) has been very successful identifying 
quasars and now there are more than 100 SDSS high-redshift quasars known (Anderson et al. 
2001). However, most of those quasars have been selected based on their color properties and 
therefore may be biased against quasars with different photometric properties like heavily 
obscured ones. Recent studies suggest that this bias may not be too problematic though 
(Richards et al. 2003). 

The other main quasar selection technique is X-ray detection as X-ray emission appears 
to be a universal characteristic of quasars at all observed redshifts (Kaspi et al. 2000). 
However, given the relative faintness of the majority of the sources, X-ray surveys need to 
rely on optical spectroscopy for most of the identifications. The lack of spatial resolution 
of past X-ray missions made these efforts laborious so far. However, the Chandra X-ray 
satellite with its superb angular resolution has made the optical identification process much 
more efficient. The nature of the X-ray emission also ensures that the optical selection biases 
due to obscuration are strongly reduced (although the X-ray selection may introduce biases 
of its own). This is especially true at high redshift where the observed X-ray photons were 
emitted at higher energies and can penetrate even considerably large amounts of obscuring 
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material. 

At high redshifts (z > 4), there are only five quasars known selected by their X-ray 
emission. Three of them were found with the ROSAT satellite (Henry et al. 1994; Zickgraf 
et al. 1997; Schneider et al. 1998). The other two were detected by Chandra (Silverman et al. 
2002; Barger et al. 2002). All of them are radio quiet except one (Zickgraf et al. 1997). There 
is another active galactic nucleus detected at z = 4.14 in the Chandra Deep Field North 
(Barger et al. 2002). However, this source is optically faint, Mb = —21.7 + 5 x Zog(/i 65 ), and 
therefore it is hard to estimate the contribution of the active nucleus itself and its putative 
host galaxy to the observed optical luminosity. We therefore exclude this object from further 
discussions. There is also an optically faint radio loud active galactic nucleus at z = 4.42 
in the Chandra Deep Field North (Brandt et al. 2001; Vignali et al. 2002c), but this object 
was selected in the radio. In total there are 56 active galactic nuclei 9 detected in X-rays 
at z > 4, 10 seven of which are radio-loud quasars, two of which are low-luminosity AGN 
(one of these is radio-loud) and twelve of which have not appeared in the literature so far. 
The remaining 36 sources are radio-quiet quasars and will be the sample we will discuss 
in this paper. Thirty-two were optically selected and four X-ray selected. Twenty-six were 
detected by Chandra, nine by ROSAT and one by XMM-Newton. These 36 sources have 
typical absorption corrected X-ray fluxes between 10~ 15 and 3 x 10 -14 erg cm -2 s _1 in the 
0.5-2.0 keV band. The absolute magnitudes of the optically selected quasars range from 
M B ~ —26 to —30 + 5 x log(h e5 ) while the X-ray selected quasars range from M B ~ —23 to 
—27 + 5 x log(h 65 ). Thus, the high redshift X-ray selected quasars sample a different range 
of the optical luminosity function than the optical selected ones. 11 

There have been several studies of high redshift optically selected quasars detected in 
X-rays (e.g., Kaspi et al. 2000; Vignali et al. 2001; Brandt et al. 2002a,b; Bechtold et al. 2002; 
Vignali et al. 2002b, c). The overall properties of the high redshift population are similar to 
the low redshift quasars (Kaspi et al. 2000; Vignali et al. 2002b). The X-ray spectral slope 
appears to be similar T ~ 2.0 at low and high redshift (e.g., Vignali et al. 2002a,b), although 
in some high redshift cases it is harder (Bechtold et al. 2002). The optical-to-X-ray spectral 
index, a ox , nonetheless, appears to be steeper (optically stronger and/or X-ray weaker) at 
high redshift than locally (Vignali et al. 2001; Bechtold et al. 2002), although this trend 



9 Here, AGN arc defined in a broad sense including quasars and active galaxies 

10 See the excellent web page http://www.astro.psu.edu/users/niel/papers/highz-xray-detected.dat main- 
tained by Niel Brandt and Christian Vignali for a list of the high redshift AGN detected in X-rays so far. 

11 Whenever we talk about optical magnitudes or luminosities we do it in a broad sense including optical 
and ultraviolet. The magnitudes in the observed frame are measured in the optical, however they correspond 
to rest frame ultraviolet absolute magnitudes or luminosities. 
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may only be a result of the possible correlation of the optical-to-X-ray spectral index with 
optical luminosity. Vignali et al. (2002d), studying a sample of SDSS Early Data Release 
quasars, find that a ox is correlated with rest-frame ultraviolet (UV) luminosity but does not 
depend significantly on redshift. They also find a correlation (with a slope different than 
unity) between X-ray and rest frame UV luminosities (see also references therein). 

Here we present the discovery of CXOCY J125304.0-090737, az = 4.179 quasar selected 
in X-rays from Chandra archival images. Optical follow-up images were taken with the 
CTIO/4- meter telescope and the confirming optical spectrum was obtained with FORS2 on 
VLT. These observations were taken as part of the Calan-Yale Deep Extragalactic Research 
(CYDER) Survey (Castander et al. 2003; Treister & Castander 2003). The CYDER survey 
is a deep optical and near-infrared imaging and spectroscopic program encompassing several 
scientific goals. One of the key aims is the characterization of the population of faint X-ray 
sources. Therefore, some of the fields of the CYDER survey were selected to overlap with 
deep Chandra pointings. 

In §2 we describe the X-ray, optical and near-infrared observations that lead to the 
discovery of CXOCY J125304. 0-090737, presenting the data reduction procedures followed. 
In §3 we discuss the observed characteristics of CXOCY J125304. 0-090737 and compare 
them to those of the high redshift population of radio quiet quasars. We also discuss the 
differences between the X-ray and optically selected high redshift QSOs. Finally we present 
our conclusions in §4. 

Throughout, we assume H Q = 65/i 65 km s _1 Mpc -1 , il a = 0.3 and = 0.7. We define 
the photon index T as the exponent giving a photon flux density X-ray spectrum f(E) oc E~ T 
in photons cm~ 2 s _1 keV -1 . 

2. Observations and Data Analysis 
2.1. X-ray Data 

The first X-ray field studied by the CYDER survey in the Southern Hemisphere was 
the Chandra pointing towards the Hickson group of galaxies HCG62 centered on (J2000) 
a = 12 ft 53 m 05.70 s 5 = -09°12'20.0" (PI: Vrtilek). This field was observed on January 25, 
2000 with the ACIS-S detector configuration for a total 49.15 ks. We retrieved this image 
from the archive and analyzed it using standard techniques with the CIAO package. 
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We remove bad columns and pixels using the guidelines specified on the ACIS Recipes, 12 , 
and we remove flaring pixels using the FLAGFLARE routine. We use the full set of standard 
event grades (0,2,3,4,6) and create two images, one from 0.5 to 2.0 keV and one from 2.0 to 
8.0 keV. We then use the WAVDETECT routine from the CIAO package to identify the point 
sources within these images, checking all wavelet scales from 1 to 16 separated by factors 
of 2. We detect 30 sources from 0.5-2.0 keV and 19 from 2.0-8.0 keV (with 15 detected in 
both bands, so that a total of 34 sources are detected) in the S3 ACIS CCD. We then use 
the PSEXTRACT script to extract spectra for these sources and provide a response matrix. 
CXOCY J125304.0-090737 was detected in both bands. Figure 2 shows a 30"x30" cutout 
region centered on CXOCY J125304.0-090737 of the total (0.5-8.0 keV) band adaptively 
smoothed X-ray image. The total number of background subtracted photons detected were 
17.1 ± 4.2 in the soft [0.5-2.0] keV band and 8.3 ± 3.0 in the hard [2.0-8.0] keV band. This 
analysis will focus on CXOCY J125304. 0-090737; the sample as a whole will be discussed in 
a future work. 

We then fit the spectrum for CXOCY J125304.0-090737 in XSPEC 11.0 with a power 
law model. We find a best fitting photon power law index of 1.57lo!33 with no absorption and 
of 1.69+035 (1°" errors for the spectral indices) when the Galactic neutral hydrogen column 
of 3 x 10 20 (as calculated using the HEASARC nH tool) is accounted for. These fits are 
done with three spectral channels from 0.5 to 8.0 keV, which all suffer from statistics too 
poor for conventional x 2 fitting and the resulting values of x 2 are close to zero, indicating 
that the bins have too few photons. Still the same results are obtained when either the C 
statistic or the x 2 statistic are used, so they can be considered as robust as any spectral fits 
based on less than 30 photons. The 0.5-8.0 keV Galactic absorption corrected observed flux 
is 4.9 ± 1.4 x 10~ 15 ergs/sec/cm 2 with either model. Figure 1 shows the observed cumulative 
energy distribution. 



2.2. Optical and Near-Infrared Imaging Data 

As part of the CYDER survey we obtained V and I images of this field with the CTIO 
4m MOSAIC-II camera. The total integration time was 80 minutes in V and 25 minutes 
in I under 1.1" seeing conditions. Images were reduced using standard techniques with the 
IRAF/MSCRED package. 13 We search for optical counterparts of our X-ray detections in 



12 See Clean the Data at http://www.astro.psu.edu/xray/acis/recipes/clean.html 

13 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
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Fig. 1. — Observed cumulative X-Ray photon energy distribution. The best fitting power-law 
model is over-plotted as a solid line, together with the 1 sigma deviations (dashed lines). 

the resulting images. The source matching is relative simple given the good spatial resolution 
of Chandra. Using the brightest nine point sources in the optical we find an overall offset 
between the X-rays and the optical (tied to the USNO reference frame) of ARA = 0.9" and 
AD EC = 1.0" (see Castander et al, in preparation for further details). There is only one 
detectable optical counterpart within the error circle of the X-ray source CXOCY J125304.0- 
090737 down to our magnitude limit (V < 25.8 and / < 23.8 at 5a). The offset between the 
X-ray and optical centroids are 1.3". For this object we performed PSF fitting photometry 
measuring V = 23.65 ± 0.05 and / = 22.52 ± 0.08. Image cut-outs of this object are shown 
in Fig 2. 

Near-infrared data for this field was obtained at LCO using the DuPont 2.5m telescope 
with the Wide Field InfraRed Camera. The total integration time was 62 minutes in the J 
band and 60 minutes in the K s band, with seeing of 0.7" and 0.6" respectively. We reduced 
the data with the IRAF DIMSUM package following standard procedures. We barely detect 
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Fig. 2.— V and X-ray images of CXOCY J125304.0-090737. Both images are cut-outs of 30" 
x 30" of the original images centered on the RA and DEC of the X-ray source. North is up 
and East to the left in both images. The V image (left) was taken at the CTIO-4m mosaic 
in a total exposure of lh 20m. The cross indicates the position of the X-ray centroid. A 
2" radius circle is also plotted. The X-ray image (right) was taken with the Chandra X-ray 
satellite with the ACIS-S chip s3, exposure time 49 ks. The X-ray image correspond to the 
full (0.5-8 keV) band and has been adaptively smoothed using the csmooth routine within 
CI AO. We apply a Gaussian kernel to the data, and use a fft convolution. The resulting 
image pixels are required to have a minimum significance of 3 a and a maximum significance 
of 5 a. 

the optical counterpart in the J band image at a magnitude of J = 21.90 ± 0.20 but we do 
not detect it in the K s image which reached a magnitude limit of K s < 19.95 (5<x). 

2.3. Optical Spectroscopy 

Follow up spectroscopy of this field was obtained at the ESO Cerro Paranal Observatory 
with the UT4/Yepun telescope using the FORS2/MXU instrument on the 14th February 
2002. Three masks were design for the HCG62 field, one of which included a slit for the 
CXOCY J 125304. 0-090737 optical counterpart. We took three exposures with this mask for 
a total integration time of 1 hour 55 minutes in ~ 0.9" seeing conditions. We used the V300 
grism giving a resolution of R ~ 520 (10.5 A) for our 1" slit. The spectrum was reduced 
using standard techniques with IRAF and calibrated in wavelength using He-Ar comparison 
lamps exposures and the night sky lines. We also obtained spectra of the spectrophotometric 
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standards Feige67 and LTT6248, which were used to provide a rough flux calibration of the 
spectrum. 

Figure 3 shows the final spectrum together with the mean spectrum of quasars at z > 4 
found in the SDSS. The typical broad emission lines of a quasar are clearly distinguishable 
and marked in the figure. We measure a mean redshift of z — 4.179 ± 0.006 using all the 
emission lines but Lyman a and NV. 
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Fig. 3.— Optical spectrum of CXOCY J125304.0-090737 obtained at UT4/Yepun VLT. The 
spectrum has been smoothed in logarithmic space with a Gaussian of 250 km/s, equal to 
the spectral resolution of the instrument for the set up used, R ~ 520 (v ~ 250 km/s). 
For comparison, we also show the error weighted average of the SDSS Early Data Release 
QSO spectra at z > 4. Each individual spectrum used in the z > 4 composite spectrum 
has a resolution of R ~ 1800 (v ~ 70 km/s). This composite spectrum has arbitrary zero 
point and scaling offsets for display purposes and it is shown above the spectrum of CXOCY 
J125304. 0-090737. The most common QSO emission lines are indicated as dotted lines. The 
emission lines in CXOCY J125304. 0-090737 are narrower than the typical SDSS spectrum. 
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3. Discussion 

Table 1 presents the observed properties of CXOCY J125304. 0-090737. In the opti- 
cal, CXOCY J125304. 0-090737 is a faint quasar with an absolute magnitude, M 1450 (i+ Z ) = 
—23.31 + 5 x log(h 65 ) and M B = —23.69 + 5 x log(h 65 ). At this absolute magnitude, one 
would expect that most of the observed flux is coming from the QSO itself with little con- 
tribution from its host galaxy, which is not resolved in our optical images. Its spectrum is 
typical of a quasar with broad emission lines (Figure 3). We measure a velocity width of 
1500 ± 300 km/s using the CIV emission line, which is lower than the one measured in the 
SDSS composite QSO spectrum (~ 2800 km/s). On the other hand, the CIV rest-frame 
equivalent width (EW = 43 ±10 A) is larger than the SDSS QSO composite spectrum (~ 24 
A) as expected from the Baldwin effect, an anti-correlation between 1450 A luminosity and 
CIV EW (Baldwin et al. 1977). The spectrum also presents a strong absorption feature be- 
tween the Lya and NV emission lines. This feature is at approximately the same wavelength 
as the [OI] 6364 A sky line and therefore its real strength is difficult to estimate. We are 
nevertheless confident that the absorption is real and not an artifact of the sky subtraction. 
We interpret this absorption line as intrinsic absorption of NV. It is approximately at 3000 
km/s blue- ward of the emission line. Although the spectrum is noisy, one may tentatively 
argue that there are also absorption lines blue-ward of CIV and OVI at approximately the 
same velocity offset. If this were correct, that would imply the existence of material that 
could obscure the QSO. Another interpretation would be that the Lya emission line is nar- 
row and that we are seeing the trough between emission lines. Although the wide red tail of 
NV would disfavor this hypothesis. The Lya flux deficit, D A) is a measure of the absorbed 
flux blue-ward of Lya (Oke & Korycansky 1982; Schneider et al. 1991). We obtain a value 
of D A ~ 0.46 for CXOCY J125304.0-090737, which is typical of quasars at this redshift 
(e.g., Madau 1995). Our measurement is only approximate due to the difficulty of fitting 
the continuum red-ward of Lya in our poor signal-to-noise spectrum, which translates into 
a large uncertainty in the extrapolation of the continuum blue-wards to compute Da- 

Using our photometry, we measure a V - 1 = 1.13 color for CXOCY J125304.0-090737. 
This value is slightly bluer than the expected V — I = 1.36 color obtained redshifting the 
SDSS composite spectrum (Vanden Berk et al. 2001) and convolving the spectrum with 
the filter responses. The similarity between these values argues against large amounts of 
intrinsic reddening (and thus optical extinction) in CXOCY J125304. 0-090737. However, 
this argument is only qualitative and does not exclude small amounts of extinction. Overall, 
considering the possible relatively weak intrinsic absorption lines in our spectrum of CXOCY 
J125304. 0-090737 and its not very red photometric properties, we conclude that there may 
be small amounts of extinction in CXOCY J125304.0-090737. 
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In X-rays, CXOCY J125304.0-090737 is also a faint QSO. The number of detected 
photon is low, 17.1 ± 4.2 in the soft [0.5-2.0] keV band and 8.3 ± 3.0 in the hard [2.0-8.0] 
keV band, and therefore fitting a power-law spectrum to the observed counts is difficult (see 
Figure 1). Nevertheless, the best photon index slope is T = 1.69^35, which is harder than 
the average radio-quiet QSO X-ray spectrum at low redshift (Yuan et al. 1998; George et 
al. 2000; Mineo et al. 2000; Reeves & Turner 2000) (although part of the difference may be 
due to selection biases) and similar to the value of the X-ray background spectrum, T ~ 1.4 
(e.g., Gendreau et al. 1995; Gilli et al. 2001). We have used this value (r = 1.69) to convert 
from observed count rate to X-ray flux. We measure an extinction corrected X-ray flux of 
4.9 ±0.5 x 10" 15 erg cm" 2 s" 1 in the observed frame [0.5-8.0] keV band. This QSO is therefore 
the second faintest X-ray selected quasar detected so far (see Figure 4) and the fourth X-ray 
faintest high redshift quasar overall. 

For spectra with low number of counts it is customary to characterize the spectral shape 
with the hardness ratio HR = (H — S)/(H + S), where S are the counts measured in a soft 
band and H are the counts in a hard band. The measured hardness ratio for CXOCY 
J125304.0-090737 is -0.35t°;i, defining the soft and hard X-ray bands as 0.5-2.0 keV and 
2.0-8.0 keV respectively. Although not many high redshift quasars have measured hardness 
ratios in the literature, the ones measured are typically softer than ours. That is their HR are 
lower and their T's are higher. Vignali et al. (2002b) measure a typical value of V ~ 2.0 ±0.2 
compared to our V ~ 1.7, emphasizing the relative hardness of CXOCY J125304. 0-090737 
compared to other optically selected high redshift QSOs. 

At redshifts z > 4, the observed Chandra band 0.5-8.0 keV corresponds to a rest frame 
band of 2.5-40 keV or higher. Any process changing the observed spectral slope is thus 
occurring at intrinsic higher energies than observed. Strong intrinsic absorption can certainly 
harden an observed X-ray spectrum. At redshift z > 4, Hydrogen column densities of 
N H > 10 23 are needed to produce noticeable spectral changes. In the case of CXOCY 
J125304. 0-090737 neither the X-ray spectrum (Fig 1) nor the optical one (Fig 3) show any 
signs of such amounts of absorption. Processes, such as Compton reflection (e.g., Guilbert & 
Rees 1988) which cause deviations from a power law spectrum at high energies may therefore 
be important for understanding the X-ray spectral properties of high redshift quasars with 
hard spectra. The most prominent such features are iron emission lines, typically seen at 6.4 
keV (from a neutral accretion disk) or 6.7 keV (from an ionized disk), iron absorption edges 
seen between 7.0 and 9.0 keV (with the energy again depending on the level of ionization in 
the disk), and a broad bump typically between 30 and 50 keV due to the Compton down- 
scattering of the hight energy photons. For typical parameters, the spectrum of a z=4.19 
quasar can be hardened by 0.10-0.15 in T with respect to the underlying power law by the 
reflection component. The effects of reflection components on the X-ray background have 



been studied in more detail in past work (e.g., Fabian et al. 1990; Zdziarski et al. 1993). 
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Fig. 4. — Soft X-ray flux (0.5-2 keV) versus ABu^i +z y Filled circles are ROSAT detections, 
empty circles are Chandra detections and the filled triangle is the only XMM-Newton de- 
tection. Squares enclosing circles represent X-ray selection, while no enclosing sign denotes 
optical selection. CXOCY J125304. 0-090737 is emphasized with another large concentric 
empty circle at A_B 1450 ( 1+2 ) = 22.95. Note that there are no optically selected QSOs detected 
in X-rays fainter than AB 1450 ( 1+Z ) ~ 20.5. The solid line represents the locus of objects with 
ot x — — 1-5 at z — 4.5 (see Fig 5). 

In order to compare the optical and X-ray emission we compute the effective optical-to- 
X-ray power law spectral slope, a ox , defined as 

_ log[/„(2 keV)/%(2500 A)] 

log[z/(2 keV)/V(2500 A)] 1 ' 

where f u is the flux density and v is the frequency of the given wavelength or energy. Neither 
flux density is measured at the given energy and wavelength and extrapolations are needed to 
compute this quantity. For CXOCY J125304.0-090737 we compute a ox = -1.35loi3 taking 
into account the error in the measured quantities and the errors in the extrapolations. 
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Fig. 5. — a x versus AB U50 ( 1+Z y Symbols are as in Fig 4: 1) small symbols denote detection: 
ROSAT, filled circle; Chandra, empty circle; XMM-Newton, filled triangle; 2) large symbols 
denote selection: X-ray, square, optical, none. CXOCY J125304. 0-090737 is highlighted with 
an extra large empty circle at ^4-Bi45o(i+z) = 22.95. The value of a ox = —1.5 is shown with 
a straight line. All optically selected QSOs have a ox values lower than -1.46. 

We have also computed a ox for all high redshift quasars (z > 4) with measured ABi45o(i+ z ) 
optical magnitudes and observed-frame 0.5-2.0 keV X-ray fluxes available in the literature. 
Figure 5 plots a ox against the optical magnitude A5i45o(i+ 2 ). The a ox values of the X-ray 
selected QSOs are shallower (X-ray stronger and/or optically fainter) than the optically 
selected ones. A Kolmogorov-Smirnov (KS) test indicates that the probability of both dis- 
tribution being drawn from the same parent distribution is less than 0.05%. The distribution 
of optical absolute magnitudes (or luminosities) is also statistically different in the X-ray and 
optically selected samples because the X-ray surveys are more sensitive but cover a much 
smaller area than the optical ones (see Figures 4 and 6). In order to understand the difference 
in the a ox values, we plot the X-ray against the rest frame UV luminosity densities in Fig- 
ure 6. We find that both observables are correlated with a functional form L 2 keV ^ L 0,55 », 

1450 A 

shown as a dashed line in Figure 6. This slope is somewhat shallower than that found by 
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Fig. 6. — X-ray .vs. optical luminosity densities. Symbols are the same as in Figures 4 
and 5. The dashed line is the best fitting linear regression including all points, which gives 
a slope of 0.55. The dotted line is a fit to the optical selected QSOs fixing the slope at 0.75 
as found by Vignali et al. (2002d). 



Vignali et al. (2002d). But taking into account that the faint UV luminosity points, driving a 
shallower slope, are X-ray selected and assuming that the distribution of X-ray luminosities 
are the same at faint UV rest frame luminosities than at bright UV luminosities, our data are 
consistent with the correlation found by Vignali et al. (2002d), thus extending their result 
to fainter luminosities at high redshift with a larger sample. 

CXOCY J125304. 0-090737 was discovered in the first Chandra X-ray Southern field 
studied within the CYDER survey. One can ask how many high redshift QSOs one would 
expect to detect given the QSO luminosity function. We use Fan et al. (2001a) QSO optical 
luminosity function (LF), which is consistent with previous estimates (e.g., Kennefick et 
al. 1995; Schmidt et al. 1995). They compute the optically selected QSO LF in the range 
3.6 < z < 5.0 and —27.75 < M U50 — 5 x log(h e5 ) < —25.75. The optical absolute magnitude 
of CXOCY J 125304.0-090737 is considerably fainter. If we extrapolate Fan et al's power law 
LF down to the magnitudes reached by our search, the expected number of quasars in the 
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ACIS-S S3 CCD with z > 4 is 0.22. The probability of finding any quasars (one or more) is 
then 20%. If instead we extrapolate Fan et al's LF to fainter luminosities using the shape 
of the optical QSO LF at lower redshifts (Boyle et al. 2000), the expected number of QSOs 
is even lower. Although, with just one detection these numbers are not significant, they are 
nevertheless indicative that the faint end slope of the QSO optical LF may be high (at least 
higher than at low redshift) or that optical surveys miss significant fractions of QSOs that 
may be obscured. In the case of CXOCY J125304. 0-090737 is worth pointing out that its 
colors are typical of low redshift quasars and therefore a mult i- wavelength optical survey 
reaching this magnitude limit should have been able to pick it up. 

The number of expected z > 4 quasars can also be computed from the X-ray luminosity 
function (XLF). Unfortunately, the XLF is poorly determined above z — 4, given the lack of 
X-ray selected QSOs at these redshifts. Compiling several X-ray surveys, Miyaji et al. (2001) 
derive an XLF up to z = 4.6, although at their highest redshift their XLF is only valid at 
the bright end. If we extrapolate their power-law parameterization of their last redshift bin 
(z = 2.3 — 4.6) to the faint end and also assume that the XLF is the same beyond z = 4.6, 
then we would expect to detect 0.72 quasars at the X-ray flux limit reached in our first 
CYDER field. If we instead extrapolate Miyaji et al's XLF with a smooth double power-law 
to the faint end similar to the one at lower redshifts and also extrapolate this function to 
higher redshift with a plausible evolutionary model (see Castander et al. in preparation 
for details), we obtain somewhat lower expectation values. The probability of finding any 
quasars would then be < 50%. We want to stress nevertheless that these numbers should be 
taken cautiously. They are only orientative to give a sense of how the faint end of the high 
redshift quasar luminosity function may behave. 

The discovery of CXOCY J125304.0-090737, together with the other X-ray selected high 
redshift QSOs (Henry et al. 1994; Zickgraf et al. 1997; Schneider et al. 1998; Silverman et 
al. 2002; Barger et al. 2002), demonstrates the potential of X-ray techniques to detect high 
redshift QSOs. At redshifts z > 4, X-ray photons observed for example at 1 keV were emitted 
at energies > 5 keV, where the X-ray photons can penetrate even considerable amounts of 
extinction and therefore X-ray selected samples should be almost extinction-bias free. Only 
sources with intrinsic column densities N H > 10 23 cm 2 should be affected. In particular, 
Chandra with its superior image quality that allows secure and economical identifications 
of optical counterparts will permit the completion of follow-up surveys, some of them which 
are already in progress, like the CYDER survey, that augur a new understanding of the 
quasar population at high redshift. For example, one can argue that one needs of the order 
of 100 high redshift quasars to constraint the high redshift LF. This could be achieved with 
40-200 Chandra pointings (ACIS-I and ACIS-S) depending on whether the actual QSOs LF 
at the faint end is as our discovery implies or lower. Although demanding, current on-going 
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X-ray follow up surveys may achieve these numbers. These discoveries will be complemented 
by deep wide-field multi-band optical surveys to provide a good understanding of the high 
redshift quasar population and thus the history of the accretion power onto black holes which 
is thought to be closely related to the galaxy formation and star formation rate history of 
the universe. 



4. Conclusions 

We present the discovery of CXOCY J125304. 0-090737, a new high redshift quasar at 
z = 4.179 found as part of the CYDER survey. CXOCY J125304.0-090737 is an optically 
faint quasar (Mb = —23.69 + 5 x log(hQ^)) with a typical QSO spectrum (Figure 3). In the 
optical, the quasar seems to be slightly obscured, although the data is also consistent with no 
extinction. In X-rays, CXOCY J125304.0-090737 is also X-ray faint (f x = 1.7 ±0.4 x 10~ 15 
ergs s -1 cm" 2 in the [0.5-2.0] keV band) with a somewhat harder spectrum (T ~ 1.7 or 
HR ~ —0.35) than typical low redshift or high redshift optically selected quasars. We 
speculate that a reflection component can slightly harden the spectrum but by no means 
is this the only mechanism. The spectrum of CXOCY J125304.0-090737 and the other X- 
ray selected high redshift quasars is X-ray strong (and/or optically weak) compared to high 
redshift optically selected quasars. Given that both samples have different optical/ultraviolet 
absolute magnitudes, this result is a manifestation of the correlation between X-ray and rest 
frame ultraviolet luminosities. Taking into account selection biases, we find a slope for this 
correlation consistent with the value obtained by Vignali et al. (2002d). We therefore extend 
their result to lower luminosities at high redshifts using a larger sample. 

The quasar space density implied by the discovery of CXOCY J 125304.0-090737 is high: 
a factor > 5 larger than reasonable extrapolations of the quasar optical luminosity function 
at high redshift. Although not statistically significant, this number seems to indicate that 
the shape of the quasar luminosity is changing with the faint end slope becoming steeper with 
redshift or that current optical surveys at high redshift are missing considerable amounts of 
quasars. The space density predicted from the X-ray luminosity function is more a guess 
than a prediction. Anyway, if on-going and future X-ray surveys corroborate the numbers 
implied by the first field of the CYDER survey, the X-ray luminosity function will not differ 
much of its value at lower redshifts, z = 2 — 3. The combination of the currently ongoing 
optical and X-ray surveys with complementary regions of the parameter space sampled will 
help to determine the quasar luminosity function at high redshift with a better understanding 
of the biases involved with both selections. The CYDER survey promises to be one of the 
contributors to a better understanding of the quasar population at high redshift. 
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Table 1. Summary of observed properties of CXOCY J125304. 0-090737 

Parameter Value 

RA (J2000) 12 fe 53 m 04.0 s 

Dec (J2000) -09°07'37" 

N H a 2.96 x 10 20 cm" 2 

z 4.179 ±0.006 

aox ®ox = _ l-35lo!i3 

V mag (Vega) 23.65 ± 0.05 

I mag (Vega) 22.52 ± 0.08 

J mag (Vega) 21.9 ± 0.2 

K s mag (Vega) > 20.5 

Mj -23.45 + 5 x log(h 65 ) 

AB 1450{1+z) 22.95 ±0.10 

Mi450(i+^) -23.0 ± 5 x log(h 65 ) 

fx (0.5-8 keV) 4.9 ± 1.3 x 10~ 15 ergs s" 1 cm" 2 

fx (0.5-2 keV) 1.7 ± 0.4 x 10~ 15 ergs s" 1 cm" 2 

fx (2.0-8 keV) 3.2 ± 1.0 x 10~ 15 ergs s" 1 cm" 2 

Hardness Ratio -0.35+[J;| 



Calculated using HEASARC tool nh 



